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Abstract: Observations of the blazars 1ES 0647+250 and 1ES 0806+524 with VERITAS are reported 
here. These objects are among the favoured candidate extragalactic sources in the very high-energy 
regime due to the presence of high-energy electrons and adequate seed photons. The presence of high- 
energy electrons is established from the location of the synrchrotron peak in the spectral energy distri- 
bution of the blazars. The presence of adequate seed photons is determined by the flux in the radio- 
through-optical wavebands. These are the key ingredients for very high-energy gamma-ray emission in 
the context of the synchrotron self-Compton model. The redshift of 1ES 0647+250 has been tentative- 
ley reported as 0.203 and the redshift of 1ES 0806+524 is 0.138, thus the detection of very high-energy 
gamma-ray emission from these objects could make significant contributions to the understanding of the 
extragalactic infrared background light. The analysis of these data relies on standard techniques in very 
high-energy gamma-ray astronomy, and the results are compared to previously reported upper limits and 
to theoretical predictions. 
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Introduction 

The blazars 1ES 0647+250 and 1ES 0806+524 are 
among several blazars observed with VERITAS 
during its commissioning phase in late 2006 and 
early 2007. The observations of 1ES 0647+250 
and 1ES 0806+524 are motivated by the search for 
very high-energy gamma-ray emission from extra- 
galactic sources. In this paper two specific models 
used to predict VHE emission from blazars will be 
discussed and previously reported measurements 
and flux upper limits of the objects compared. The 
VERITAS instrument is briefly described and the 
results from the analysis of the VERITAS data are 
reported. 

Predictive Models 



Compton components of their spectra. The sec- 
ond modification extends the radio range below 
10 41 ' 2 erg s _1 by modeling Compton scattering in 
the Klein-Nishina regime. This modification also 
uses a different width for the parabola representing 
the Compton peak, which is reduced with respect 
to the synchrotron peak. 

In the Costamante [4] model, a single zone SSC 
fit to multiwavelength data is used to predict fluxes 
in the TeV regime. The model emphasises the re- 
quirement of both high-energy electrons and suf- 
ficient seed photons to produce very high-energy 
gamma rays. In the study, a large sample of 
BL Lacs were examined and the radio/optical flux 
and synchrotron peak frequency fit using the SSC 
model. The flux predictions are summarised in Ta- 
ble 1. 



In the modified-Fossati model [4, 7], the peak fre- 
quency of the synchrotron spectrum and the rela- 
tive importance of the inverse-Compton power are 
determined by the radio luminosity. The first mod- 
ification assumes that objects of low power have 
equal luminosities in the synchrotron and self- 



Targets 

The blazar 1ES 0647+250 was discovered in the 
MHT-Green Bank survey at 5 GHz using the 
NRAO 91-m transit telescope. X-ray emission was 
discovered in the Einstein Slew Survey with the 
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Object 


Costamante 


Modified-Fossati 


1ES 0647+250 


0.59 


0.24 


1ES 0806+524 


1.36 





Table 1 : Flux predictions according to the Costa- 
mante and modified-Fossati models. Fluxes are in 
units of 10~ n ever 2 s _1 above 0.3 TeV. 



synchrotron peak falling just below 10 keV. A red- 
shift of 0.203 has been tentatively reported. Previ- 
ous observations of this object were reported by 
HEGRA [1] with a 99 % confidence flux upper 
limit of J 7 e>o.78ToV < 33.5 x 10~ n cm~ 2 s _1 
from 4.1 hours of observations, where a spec- 
tral index of -2.5 was assumed. To compare 
this to the Costamante and modified-Fossati mod- 
els, this can be extrapolated to J 7 e>o.3ToV < 
126.4 x 10 _11 cm _2 s _1 assuming a constant spec- 
tral index of -2.5 down to 0.3 TeV. This is well 
above both the predictions of the Costamante and 
modified-Fossati models and does not constrain ei- 
ther model. 

The blazar 1ES 0806+524 was discovered using 
the NRAO Green Bank 91-m telescope at 4.85 
GHz with X-ray emission reported by the Einstein 
Slew Survey. The galaxy has a measured red- 
shift of 0.138. Whipple [9, 6] reported flux upper 
limits of J 7 e>o.3ToV < 1.37 x 10~ n cm~ 2 s~\ 
^E>o.3TcV < 16.80 x lO" 11 cm" 2 s" 1 and 
^ r E>o.3TcV < 1.47 x 10 _11 cm _2 s _1 from differ- 
ent observing seasons. HEGRA reported an upper 
limit of J r E> i.o9 TeV < 42.5 x 10~ n cm -2 s _1 in 
one hour of observations which can be extrapolated 
to^E>o.3ToV < 255.4 x 10 _11 cm _2 s _1 (assum- 
ing a spectral index of -2.5) to compare with the 
Whipple result and Costamante prediction. Nei- 
ther the Whipple nor the HEGRA results constrain 
the Costamante model. 



Observations 

All observations were taken in WOBBLE mode. In 
this mode, the target is offset from the center of 
the field of view by ±0.3° or ±0.5° in declina- 
tion (or right ascension). One of the trade-offs 
with the WOBBLE mode is that although more time 
can be spent with the observational target in the 



field of view, the target is not at the center of the 
field of view where the sensitivity is highest. Thus 
the WOBBLE offset must be carefully chosen using 
measurements or Monte Carlo simulations of the 
detector response to maximise sensitivity. All ob- 
servations were taken during the VERITAS com- 
missioning phase, where initially two telescopes 
and later three telescopes were complete. Af- 
ter selecting data to remove runs suffering from 
bad weather or technical problems, a total of 17.3 
hours on 1ES 0647+250 and 34.7 hours on 1ES 
0806+524 were available. 



Data Analysis 

The data have been analysed using independent 
analysis packages (see [3] for details on the anal- 
ysis package). All of these analyses yield con- 
sistent results. Standard data analysis techniques 
[5] for ground-based gamma-ray astronomy were 
used, and are briefly described here. The analysis 
was optimised using Crab Nebula [2] data from the 
same period. 

Shower images in the focal plane are gain cor- 
rected and cleaned before being parameterised us- 
ing a moment analysis. For each shower, the fo- 
cal plane images are parameterised using Mean- 
Scaled Width and Length (MSW/MSL)[5]. These 
are calculated using a large data set of Monte Carlo 
simulations at discrete zenith (6) angles (with in- 
terpolation in cos 9). Cuts on MSL and MSW are 
designed to reject most of the background while 
retaining a large portion of the signal. The cuts for 
this analysis were optimised on a data set of Crab 
Nebula observations. 

Background estimation is performed using the 
reflected-region background model. In this 
scheme, an integration region is placed around the 
putative source position, with identical background 
integration regions distributed around the field of 
view. The number of events in the integration re- 
gion is termed On, the number of events in the 
background region is termed Off and the ratio of 
the integration areas is termed a. 

The analysis results are summarised in Table 2, 
with the data from both objects broken down in 
terms of the number of participating telescopes. 
The statistical excess is calculated using equation 
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17 from [10]. No excess above 5 a is uncovered 
in any of the data. A distribution of 9 2 for the 
signal and background regions for both sources is 
shown in Figures 1 and 2. The 9 2 distribution 
shows the squared angular distance between the 
putative source position and reconstructed shower 
source. The red line marks the 9 2 distribution rel- 
ative to the putative source position and the blue 
points mark the averaged 9 2 distribution relative to 
the background positions. The signal integration 
region is to the left of the vertical line. 

In the absence of a clear signal greater than 5 a, 
upper limits for both objects can be derived. This 
is done in order to compare the predictions of the 
Costamante and modified-Fossati models and to 
previous reports. The upper limit is found from 
the probability density function of the number of 
counts from the putative source [8]. 
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where p is the confidence level, C is the number of 
excess counts, a is given by 



a = y/ on + a 2 Off 
and the function / is given by 



I(z) 



1 



/2tt 



e~ x2 ' 2 dx 



(2) 



(3) 



The energy threshold (peak energy response after 
gamma-ray selection cuts) for each data set is cal- 
culated using a data base of Monte Carlo gamma- 
ray simulations at multiple zenith angles. The en- 
ergy threshold is found by interpolation in cos 
using the mean elevation of each data set. 

After calculating the 99 % upper limit on the num- 
ber of counts. It is then expressed in terms of Crab 
Nebula units, and converted into a flux upper limit 
by scaling relative to the Crab Nebula flux. These 
results are shown in Table 2. The differential flux 
upper limits are shown in Figure 3. 

Discussion 

The upper limits calculated for 1ES 0647+250 are 
significantly better than those reported previously. 




Figure 1: 6 2 distribution for 1ES 0647+250 - the 
signal integration region is to the left of the verti- 
cal line. The background counts are given by the 
crosses. 
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Figure 2: 6 2 distribution for 1ES 0806+524 - the 
signal integration region is to the left of the verti- 
cal line. The background counts are given by the 
crosses. 
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Figure 3: Spectral energy distributions of 1ES 
0647+250 and 1ES 0806+524 from [4], with the 
differential flux upper limits from this work over- 
laid in blue. The vertical axis displays log(vFv) in 
units of erg cm~ 2 s _1 . 
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Target 


Data Set 


Exposure (hrs) 


Excess 


a 


E TcV 


•'EX). 3 ToV 


1ES 0647+250 


2-Tel 


16 


38 


1.5 


0.275 


0.58 


1ES 0647+250 


3-Tel 


1.3 


13 


2.1 


0.225 


1.3 


1ES 0647+250 


All 


17.3 


51 


2.0 


0.275 


0.62 


1ES 0806+524 


2-Tel 


8.5 


18 


1.3 


0.345 


0.42 


1ES 0806+524 


3-Tel 


26.2 


67 


2.1 


0.36 


0.15 


1ES 0806+524 


All 


34.7 


85 


2.5 


0.36 


0.13 



Table 2: Summary of the VERITAS observations of 1ES 0647+250 and 1ES 0806+524 from November 
2006 to March 2007. The data are categorised in terms of the array configuration. The energy threshold, 
E^ v , is calculated for the average zenith angle of the dataset in TeV The Excess is On — a(9j/f and a is the 
statistical significance as calculated using Equation 17 from [10]. ^Fe>o.3 ToV i s me ^ ux u PP er limit above 
0.3 TeV in units of 10" 11 cm~ 2 s" 1 . 



The Costamante model is surpassed, however with 
an extremely low prediction of ^ 7 e>o.3ToV = 
0.24 x 10~ n cm~ 2 s~\ the modified-Fossati pre- 
diction is not surpassed. The upper limits cal- 
culated for 1ES 0806+524 are also significantly 
better than those reported before, and surpass the 
Costamante model. 

Although predictions for the Costamante model on 
both objects are surpassed, the implications are 
not necessarily severe for SSC modeling, espe- 
cially given the uncertainties in scaling of flux lim- 
its when the spectrum is unknown. The prediction 
of TeV fluxes is extremely sensitive to small vari- 
ations in model parameters such as radio and X- 
ray flux. Also, the flux predictions were built from 
multiwavelength data on sources that are known to 
exhibit variability. Thus the flux at the lower en- 
ergy regions of the data sets may have been dif- 
ferent during these VERITAS observations, than 
it was when the multiwavelength data from which 
the model was built. 

Given the lack of detectable quiescent emission 
with a modestly deep exposure, these objects are 
unlikely to be the subject of future deep obser- 
vations with an instrument of this class. How- 
ever, given that optical and x-ray flaring have been 
linked to detectable increases in emission in the 
TeV regime, future campaigns may revolve around 
target-of-opportunity triggers from the multiwave- 
length community. 



Conclusions 

No evidence for emission is found from observa- 
tions of 1ES 0647+250 and 1ES 0806+524. Con- 
straining upper limits were placed on the emission 
from 1ES 0806+524 in the TeV regime. The VER- 
ITAS blazar key science project promises exciting 
results in the future as the array reaches maturity. 
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